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Abstract 

In this thesis, the following type Tikhonov regularization problem will be system- 
atically studied: 

(u t , v t ) := argmin {\\v\\ x + t\\u\\ Y }, 

u+v=f 

where Y is a smooth space such as a BV space or a Sobolev space and X is the space 
in which we measure distortion. Examples of the above problem occur in denoising in 
image processing, in numerically treating inverse problems, and in the sparse recovery 
problem of compressed sensing. It is also at the heart of interpolation of linear operators 
by the real method of interpolation. We shall characterize the minimizing pair (ut,vt) 
for (X, Y) := (^(O), BV(fi)) as a primary example and generalize Yves Meyer's result 
in [11] and Antonin Chambolle's result in [6j. After that, the following multiscale 
decomposition scheme will be studied: 

Uk+i ■= argmin {- 1|/ - + t k \u - Ufc|Bv}, 
«eBV(Q)nL 2 (tt) 1 

where uq = and Q is a bounded Lipschitz domain in M. d . This method was intro- 
duced by Eitan Tadmor et al. and we shall improve the L2 convergence result in |16j . 
Other pairs such as (X, Y) := {L p ,W l {L T )) and (X,Y) := (£2,^) will also be men- 
tioned. In the end, the numerical implementation for (X, Y) := (£2(^)5 BV(fJ)) and 
the corresponding convergence results will be given. 

1 Introduction: The Importance of Research 

Many problems in optimization and applied mathematics center on decomposing a given 
function / into a sum of two functions with prescribed properties. Typically, one of these 
functions is called a good function u and represents the properties of / we wish to maintain 
while the second part v represents error/distortion or noise in the stochastic setting. Exam- 
ples occur in denoising in image processing, in numerically treating inverse problems, and 
in the sparse recovery problem of compressed sensing. The general problem of decomposing 
a function as a sum of two functions is also at the heart of interpolation of linear operators 
by the real method of interpolation. My research explores the mathematics behind such 
decompositions and their numerical implementation. 



1 



One can formulate the decomposition problem for any pair of Banach spaces X and Y 
with Y the space of good functions and X the space in which we measure distortion. Given 
a real number t > 0, we consider the minimization problem: 

K(f,t):= inf {\\ v \\ x + t\\u\\ Y }. (1.1) 

f=u+v 

K(f,t) is called the K-functional for the pair (X, Y). The pair (u t ,v t ) which minimizes 
K(f, t) is the Tikhonov regularization pair: 

(u t ,v t ) := argmin {|MU + *|M|y}. (1-2) 

u+v=f 

One can usually prove (by compactness argument and strict convexity) that there exists a 
unique solution (u t ,v t ) for problem (11.21) . As we vary t, we obtain different decompositions. 
These decompositions describe how / sits relative to X and Y. There are many variants of 
(ll.2p that are commonly used. For example, the norm of Y can be replaced by a semi-norm 
or a quasi-norm and sometimes the norm with respect to X is raised to a power. 

While the above formulation can be defined for any pair (X, Y) and any / in X + Y , 
in applications we are interested in specific pairs. One common setting, and the first one of 
interest to me, is when X = L p and Y is a smooth space such as a Sobolev space or a BV 
space. This particular case appears in many problems of image processing, optimization, 
compression, and encoding. We shall study various questions associated to such decomposi- 
tions. 

The main problems to be investigated in this thesis are: 

(i) Given / and t > 0, characterize the minimizing pair (u t ,v t ). 

(ii) Find an analytic expression for K(f, t) in terms of classical quantities and thereby 
characterize the interpolation spaces for a given pair (X, Y). 

(iii) Multiscale decompositions corresponding to the pair (X, Y) that can be derived from 
the characterization of the minimizing pair. 

(iv) Numerical methods for computing this decomposition or something close to it. 

The structure of this thesis is as following: 
Chapter 1: Introduction: The Importance of Research. 
Chapter 2: Basic properties of BV(fi) and Hausdorff measure. 

Chapter 3: Decomposition for the pair (L 2 (fi), BV(fi)), where Q is a bounded Lipschitz 
domain in M d . In this section, we first characterize the minimizing pair (u t ,v t ) by 
studying the Euler-Lagrange equation associated to (11.21) which includes giving an 
appropriate setting for the boundary condition. We generalize Yves Meyer's result 
in [TT] and Antonin Chambolle's result in [6] on the properties of (ut,Vt). Then the 
expression of K(f, t) follows as a simple consequence. In addition, we propose simpler 
proofs about characterizing the sub differential of BV semi-norm which were first proved 
in ®. 
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Chapter 4: Multiscale decompositions corresponding to the pair (Z/2,BV). In this section, 
we study the scheme introduced by Eitan Tadmor et al. under the general framework 
of Inverse Scale Space Methods and improve the L2 convergence result in [16] . 

Chapter 5: Decomposition for (X,Y) := (L^W 1 ^)) with 1/r := l/p + l/d. 

Chapter 6: Decomposition for (X,Y) := (£2,i p ) with 1 < p < 00. 

Chapter 7: Numerical implementation for (X,Y) := (L2(fi), BV(fi)) and the correspond- 
ing convergence results. 

2 The Space BV(fi) and Hausdorff Measure 

In this section, we will introduce some basic facts about the space BV(fi), where Q is a 
bounded Lipschitz domain. First, we need to give a definition of the Lipschitz domain. 

Definition 2.1 (Lipschitz Domain). An open set Q C M. d is a Lipschitz domain if for any 
xq G dQ there exists r > and a Lipschitz function $ : — > R such that - upon relabeling 
and reorienting the coordinates axis - we have 

Q n B(x , r) = {x G B(x , r) : x d > . . . , 

In the following text, without specifically mentioned, we will assume Q is a bounded 
Lipschitz domain in R d and use the following notations: 

• Lebesgue measure of Q. 

• x := (xi,x 2 , • • . ,x d ). 

V 9xi ' 8x2 ' ' ' ' ' Sa;^ / ' 

• D a u := d%ld%* ■ ■ ■ <9^w, where a := a 2 , . . . , «d) and «j G N U {0} for 1 < z < d. 

Before introducing the space BV(fi), we need to give definitions of weak derivative and 
measure. 

Definition 2.2 (Weak Derivative). Let u G Li(fl). For a given multi-index a, a function 
u G is called the a th weak derivative of u if 

/ v<\)dx= (-1) |Q| / uD a (f) dx 
Jn Jn 

for all G Co°(fi). Here |a| := «i + a 2 + • • • + a^. 
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Definition 2.3 (Measure). The a th weak derivative of u is called a measure if there exists 
a regular Borel (signed) measure /ionfi such that 

f <j>dn= (-1) H I uD a <p dx 
Jn Jn 

for all G C£°(Q). In addition |D Q -u|(fi) denotes the total variation of the measure \i. 
Definition 2.4. 

BV(fi) := Li(fi) n {u : -D a w is a measure, |.D Q u|(fi) < oo, |or| = 1}. 
In addition, the BV semi- norm |m|bv can be defined as: 

| M | BV = j \Du\ : = sup { / wdiv(0) dx : <p G C™{Q;R d ), \<f>(x)\ < 1 for x G Q}. 
Jn Jn 

Remark. It is easy to see, for u G W 1 (Li(D,)) , |m|bv = Jo |V«| cfcc. 

Theorem 2.5 (Coarea Formula). Lei w G BV(fi) and define E t := {x £ Q : u(x) < t}. 
Then 



r- poo r 

/ \Du\ = dt \DX Et \ 
J n J — oo J n 



Proof. See Theorem 1.23 in [Hj. □ 

Definition 2.6 (Hausdorff Measure). For set C W 1 , < k < oo and < 5 < oo, we 
define 

oo oo 

H s k {E) := u k 2~ k inf {^(diam Sj) k : E C (J S,-, diam ^ < 5} 

i=i i=i 

and 

H\E) :=KmflJ(J5) = supfl][(JE0, 

<5-».0 5 

where Co>fc = r(|) fc /r(| + 1), > 0, is the Lebesgue measure of the unit ball in M. k . 
H is called the A;- dimensional Hausdorff measure. 

Example 2.7. Suppose E dVt has C 2 boundary and consider Xe the characteristic function 
of E, then 

\x e \bv = H d -\dEnn). 

Proof. See Example 1.4 in [8]. □ 

For a bounded Lipschitz domain Q, the outward unit normal vector u(x) := (z/ 1; u 2 , Vd) 
is defined H d ~ l -&.e. on <9f2. Then we have the following generalized Gauss-Green theorem: 

/ div(F) dx= [ F(x) ■ u(x) dH^ 1 
Jn Jan 

whenever F G C 1 (fi;IR d ). For detailed exposition, please refer to [T7] . 
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3 Decomposition for (X, Y) = (L 2 (0), BV(fi)) 



While there are many settings and potential decompositions that we shall discuss, a partic- 
ular problem which is of high interest and is a primary example of the goal of my research is 
the problem of decomposing a function / G /^(fi), where Q is a bounded Lipschitz domain 
in IR d , into a sum of an L2 function and a BV function. For any prescribed t > 0, we define 
the pair (ut,Vt) as the solution of the following minimization problem: 

(u t , v t ) := argmin {- ||u||? + i|u| B v}- (3-1) 



u+v=f 
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If we define T(w) := |||/ — + t|w|BV, then problem (13. ip is equivalent to the following 
problem: 

u t := argmin T(u). (3.2) 

ueBV(n)nL 2 (n) 

Problem (13 .2p is closely related to the following constrained minimization problem: 

min J(u) subject to If dx = u dx and II f — u\\l = er, (3.3) 
ueBV(n)nL 2 (n) v ; J n J J n " " 2 V ; 

where J(u) := |m|bv- It is widely used in image denoising where it is called Rudin-Osher- 
Fatemi model for Q C M 2 (see [15J). If / is a given noisy image, then u t captures the main 
features of / and v t contains the oscillatory patterns of texture or the inherent noise in the 
image. In the stochastic setting, a central question is what is the best choice of t. 

Before introducing our main results, we need to spend a few words on the rigorous 
definition of the solution of (13. 2p . 

To derive the Euler-Lagrange equation associated to problem (13. 2p . we first consider a 
special case that u G C 1 (fi) and dVL is C 1 . Set T{u) := — + tJ(u) with J(u) : = 
J n |Vm| dx. Consider the following minimization problem: 

mm T(u), 

ues 

where 

— Oil 

S:={ue C\tt) : — = for x G dfl and Vw 7^ for all x e Q} 
is the admissible set. 

Notice V(|x|) = o for x 7^ 0. We can thus calculate the Gateaux derivative of J{u) for 
u E S. Given u G 5, for any /i G C 1 (fi) with |£ = on when e small enough, we have 



5 



u + eh G S. Hence: 



5J(u;h) = lim-{J(w + eh) - J(u)\ 

= lim- / {\V(u + eh)\ - |Vw|} dx 
e ./o 



lim- [ {V(\x\)\ x=Vu -eVh + 0(e 2 )} dx 
e J n 



• v h dx 



n Vu 



-div( -r— )/?- dx + f div( r— - fe) 



/" ,. / ^ M m , /" h du 1 
I — divi-— — r)h dx + / r— , — — ds 



-divf r— )/i dx, 
\vu\ 

where I 5 * := V« ■ za 

Of 

The reason why we choose Neumann boundary condition is to make j n u dx = J n f dx, 
which means the error/distortion v = f — u has mean value zero. As we shall show below, 
j n u dx = j n f dx will automatically be satisfied when u is a minimizer for problem (13.21) . 

The necessary condition for u to be a minimizer is: 5T(u; h) = which means 

t5J(u; h) — / (/ — u)h dx = 



for any h G C 1 (fi) with ^ = 0. Hence, we can informally write the Euler-Lagrange equation 
associated to problem (13. 2p as: 

Ti - *div(^r) = / infi 



[V«[> 

I s = o on an 

01/ 



Now we come back to the more general case which u G BV(fi) H L 2 (fi) with 

J(u) := |m|bv = sup { / wdiv(0) dx : <j> G V^}, 

in 

where ^ := {0 G C^O; M d ) : \(j)(x)\ < 1 for x G fi}. 

We can extend the domain of to L 2 (fl) in the following way: 



J(u) := 

In this way, we can also define 



Mbv u G BV(fi) n L 2 (Q) 

+oo neL 2 (fi)\(BV(fi)ni 2 (n)) 



r(«):=5ll/-«lli a + ^(«) 



on ^(fi). In the following text, without specific mention, we will assume J(u) and T(u) 
defined on the whole space of L 2 (0) as above. It is easy to check that J(u) and T{u) defined 
in this way are proper convex functionals on L 2 (fi). 
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Lemma 3.1. J(u) is weakly lower semi- continuous with respect to L p (l < p < oo) topol- 
ogy, i.e. if u n — ^ u weakly in L p , we have: J(u) < lim inf n _>.oo J(un) ■ ^ n addition, if 
liminfn^oo J{u n ) < oo, we have u G BV(ft). 

Proof. For any G V, where V := {<p G Cg(ft; M d ) : |0(a;)| < 1 for x G ft}, we have 



/ wdiv(0) c?x = lim / u n div((p) dx 
Jn MO % 

= lim inf / w n div(0) dx 
n ^°° Jn 



< lim inf J(u n ). 

n— >oo 

Take over the set V, we get J(u) < lim inf n ^oo J{u n ). □ 

Now we give the existence and uniqueness result for problem (13. 2 j) without invoking the 
associated Euler-Lagrange equation. 

Theorem 3.2. Fort > 0, there exists a unique minimizer u t for problem A3.2fy . In addition, 
u t solves (TX3p for cr = ||/-w t |U a ^ II/-/IU2 an d f Q u t dx = f n f dx, where f = p f n f dx. 

Proof. Let u n be a minimizing sequence for T(u). Since ||w||l 2 < II/IU2 + V^( u )i {ll^nll^} 
are bounded. Since L p is reflexive when 1 < p < 00, there exists a subsequence {u nj } such 
that M nj — Ut weakly in L2. By the weakly lower semi-continuity of J(u) as in lemma I3TT1 
we have: 

tJ(u t ) + 7:\\f ~ u t \\l 2 < lim inf {tJ(u n ) + - \\f - u n \\l 2 } = min 

Hence ut solves (13 ,2ft and the uniqueness of the minimizer follows immediately from strict 
convexity of T{u). 

Suppose o = \\f -u t \\ L2 > \\f - /|| X2 , then 

tJ(f) + l\\f-f\\l<\\\f- ut\\l 2 < Ufa) + \\\f- WtllL, 

which is contradictory with the definition of the minimizer. So we have a = \\f — Ut\\L 2 < 
\\f-fh 2 . 

Similarly, if f n u t dx 7^ f n f dx, let c = p — u t ) dx. Then we have ||/ — {u t + 
c )l|i2 < 11/ ~~ n t 1 1 i 2 while J(w t + c) = J{u t ), which is contradictory with the definition of the 
minimizer. □ 

Remark. There is an alternative way to get the existence proof by compactness argument 
for BV(ft) (See [T]). 

To characterize the minimizer of problem (13. 21) . we have to come back to the PDE ap- 
proach. The associated Euler-Lagrange equation can be formally written as: 

u - tdiv(^) = / in ft . , 

|K = on 9ft 1 J 

01/ 

To understand this equation correctly, we first need to give a rigorous definition of the 
boundary condition and the nonlinear operator — div(p^). 
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3.1 Definition of Neumann Boundary Condition 

Throughout this section we frequently make use of results shown by Anzellotti in [3]. To 
define the Neumann boundary condition in the sense of trace, we shall consider the following 
spaces: 

• BV(fi) p := BV(fi) nLp(ft). 

• X(Q) g = {z G ^^(fijIR^) : div(z) G L q (Q)}, where div(z) is defined in the sense of 
distribution: j n div(z)(p dx = — f n z • V0 dx for any G Cq°(CI). Here H^Hl^ : = 

II \/X/i=l ^lUoo- 

Set 

{OO J9 = 1 
p3i 1 <P < oo 
1 p = OO 

If z G and u> G BV(fi) p , we can define the functional (z,Dw) : Co°(fi) i— > K by the 

formula: 

((z, Dw),(j)} :— — / w;0div(z) da; — / • V0 (ix. 



Theorem 3.3. TTie functional (z,Dw) defined above for z G X(f2) 9 and w G BV(f2) p zs a 
Radon measure on Q, and 

z, Dw) = z ■ Vw dx 



Jn 

form G W /1 (L 1 (f2))nL p (fi). In addition, we have \ f B (z, Dw)\ < j B \(z, Dw)\ < \\z\\ Loo § B \Dw\ 
for any Borel set B dVt. 

Proof. Given w G BV(f2) p , we can find a sequence {w n } C D BV(f2) p (see [8]) such 

that: 

w n — > w in and limsup / |Z?io„| < / \Dw\ for any open set Acfl. 

n-i-oo JAno J Ann 

Take any G C^(A) and consider an open set V such that supp(0) C V CC A, then we 
have 



((z, Dw n ), 0) — — w n div(cj)z) dx = / <pz ■ Vw n ckc. 
Jv Jv 

So 

|((z, £)«;„), 0)| < ||0||zoq(v) lkl|x„,(V) / |Vw„| dx. 

Jv 

Taking the limit for n — > oo, we get 

|<(z, Dw), 0}| < 11011^(^11^11^^) /_ |2>w| < IWUJN^ / \Dw\ forany0GC o co (A). 

Jvnn J A 

So (z,Dw) is a Radon measure and we have | j B (z, Dw)\ < f B \(z,Dw)\ < H^Hz,^ J B \Dw\ 
for any Borel set Bcfl. 

Since ((z, Dw), 0) = — f n u>div(0z) dx = f u 0z ■ Vw dx for any w G W /1 (L 1 (fi)) fl 
and G C °°(fi), we get L(z, Dw) = Lz-Vwdx for iu G W^L^tt)) n L p (fi). □ 



Theorem 3.4. Let z G X(Q) q , w G BV(Q) P and (z,Dw) be defined as in Theorem \ 3.3[ 
Then there exists a sequence {u> n }^l C C°°(Q) fl BV(Q) P such that 

w n — > w in L p {Vi) and / (z, Dw n ) — > / (z,Dw). 

Jn Jn 

Proof. For any e > 0, take an open set A C ft such that 

\Dw\ < e, 

Q\A 

and let g G C^°(Q) be such that < g(x) < 1 in Q and = 1 in A. We can find a 
sequence {w n } C D BV(fi) p (see [8]) such that: 

w n — > w in L p (Vl) and limsup / |Vw n | dx < \Dw\. 

n^co Jn\A Jn\A 



Then 



where 



[z, Dw n ) — / (z, Dw)\ < \((z,Dw n ) } g}-((z,Dw) } g}\ 
Jn 

+ [ \(z,Dw n )\(l-g) + [ \(z,Dw)\(l-g), 
Jn Jn 

lim ((z,Dw n ),g) = ((z,Dw),g), 

n— >oo 

limsup / \(z, Dw n ) |(1 -g) < \\ z\\ Loo limsup / |Diw n | < e||z|| Loo , 

n— >>oo Jq n— >>oo J n\A 

[ \(z,Dw)\(l-g)< f \(z,Dw)\<e\\z\\ Loo . 
Jn Jn\A 

So the theorem is proved, as e is arbitrary. □ 
Theorem 3.5. There exists a linear operator^ : X(Q) q h-> Loo(<9f2) swc/i t/icrf 

ll7(^)IU»(an) < IklUoc(n)- 

2. j(z){x) = z(x)-v{x) H d - x -a.e. on dtt for z e C\tt;R d ) . 

3. (z,w)qq := f n wdiv(z) dx + f Q (z,Dw) = f gci j(z)tr(w) dH d ~ l for any w G BV(Q) P , 
where tr(w) G L>i(dQ) is the trace of w on dfl. 

Proof. Let (z, w)sn '■= f n wdiv(z) dx+ f n (z, Dw), first we want to show (z, wi)sn = (z, w 2 )qci 
for any tr{w\) = tr{w 2 ) and W\,w 2 G BV(Q) P . We can find a sequence of functions {g n } C 
C^(fl) such that g n — > Wi — w 2 in L p and f n (z, Dg n ) — > f n (z, D{wi — w 2 )). Then we have: 



(z,Wi - w 2 )dn = / (wx - w 2 )div(z) dx + (z,D(w 1 -w 2 )) 
Jn Jn 

= lim { / g n dw(z) dx + (z,Dg n )} = 0. 
n ^°° Jn Jn 



So (z,iui)an = (z, u> 2 )an- 

Now we want to show \(z,w)qq\ < H^Hz^ J dn \tr(w)\ dH d ~ l . For bounded Lipschitz 
domain Q, given any u G L\(dQ) and e > 0, we can find a function u> G such 
that 

tr(w) — u, I | Vu>| dx < \u\ dH d ~ l + e, w(x) = for x G 
Jn Jan 

where f2 e := {x G : dist(x,<9fi) > e}. Then for any tr(w) G Li(Q), we can find w G 
VF 1 (L 1 (n)) such that tr(w) = tr(w) with the above properties. So 



(s,w)an| = \{z,w) 



an\ 



< | / wdiv(z) dx\ + II^Hl^ / \Dw\ 
<n Jn 



< 



/ wdiv(» dx\ + \\z\\ Loo { / \tr(w)\ dH^ 1 + e}. 
Jn\n e Jan 



>n\n 

Since lim e ^ Jq\ q wdiv(z) dx = 0, let e goes to 0, we get 

\(z,w)en\ < IklU- / 1^0)1 (3-5) 



Jan 

Now given a fixed 2; G X(Q) q , we can define the linear functional F z : L\(dQ) 1— > M by 

:= (z,w) d n, 

where tr(w) = w. From (13.51) . we know < ||^IUoo(n)ll u lki(9n)- By Ri eze Representa- 

tion theorem, there exists j(z) G Loo(<9fi) such that 



F z {u) = [ 7 (z)u dtf"- 1 . 
Jan 



So (z,w)an = / 9n 7(z)fr(w) 1 and 1(7(2) || Loo (an) < IMIwn)- 

When ^ G C 1 ^;^), (z,w) dQ = J n div(wz) dx = J dn tr(w)z ■ v dE d ~ x . So 7(2) = 
z ■ v H d ~ l -a,.e. on dVt. □ 

Thus the function 7(2) is a weakly defined trace on dfl of the normal component of z, we 
shall denote 7(2) by [z, z/]. In this way, the Neumann boundary condition can be expressed 
as [z, u] = H d ~ 1 -a.e. on d£l. 



3.2 Definition of the Operator — div(p^r) 

Let A : X 1— >■ 2 X be a multivalued mapping defined on a Banach space X, i.e., A assigns to 
each point u G X a subset of X*, where X* is the dual space of X. In this paper we will 
simply call such mapping an operator. 

1. The set D(A) := {u G X : Au 7^ 0} is called the effective domain of A. When 
D(A) 7^ 0, we say A is proper. 
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2. The set R(A) := {J ueX ^ u ^ s cane d the range of A. 

3. The set G(A) := G X x F : u e D(A),v G Au} is called the graph of A. In 
this paper, we briefly write (u, v) G A instead of (u, v) G G(A) and we will identify an 
operator A with its graph G(A). 

4. An operator A is called a monotone operator, if (vi—v 2 ,Ui—u 2 ) > for («]_, i^), (w 2 , t> 2 ) G 
A. 

5. A monotone operator A is called a maximal monotone operator, if for any monotone 
operator B that A C -B, we have A = B. 

Lemma 3.6. Let X be a Hilbert space and A : -D(A) C X i— >• 2 X . If A is a monotone 
operator, then for any A > ; (J + AA) -1 : X i— >• -D(A) zs nonexpansive, i.e. for V\,v 2 G 
D((/ + AA)- 1 ) ; we W \{I + \AY\v x )-{I + \A)-\vl)\x< IK-^IU- 

Proof. Let fi G (/ + AA)(wi) and v 2 G (/ + AA)(w 2 ). Then ifi = — ^i) G A(/Ui) and 
w 2 = \( v 2 — u 2 ) G A{u 2 ). By the fact (wi — w 2 , U\ — u 2 ) > 0, we have: 

\\vi-v 2 \\ 2 x = ((J + AA)( Ml )-(J + AA)(« 2 ),(J + AA)( Wl )-(/ + AA)(« 2 )) 
= ||«i - u 2 \\ 2 x + A 2 ||wi - u^Hx + 2A(wi - w 2 , u\ - u 2 ) 
> \Wi-u 2 \\ 2 x . 

□ 

Now we introduce the following operator A on L 2 (Q): 

Definition 3.7. v G A(u) means: 

3 z G X(f2) 2 with H^Hz^ < l,v — — div(z) in T>'(Q) such that 
j n (4> - u)v dx = j n (z, D<f>) - |u|bv for any <p G BV(fi) n L 2 (fi). 

Let's recall the set S that we used for deriving Euler-Lagrange Equation: 

S :={ue C\Q) : — = for x G dQ and Vm 7^ for all x G fl}. 

For u E S, A(u) = — div(j^r). Hence the operator „4 can be viewed as generalization of 
-div(^j). Formally, we can write A{u) = -div(^) for u G BV(fi) n L 2 (fi). 

To associate the operator ^4 with our minimization problem (13. 2p . we need to introduce 
the concept about subdifferential. Let L : X 1— > [—00, +00] be a functional on a real Banach 
space X. The functional u* in X* is called a subgradient of L at the point w if and only if 
L(u) 7^ ±00 and 

L(w) > L(u) + (u*,w — u) for any w G X. 

For each u e X, the set: 

dL{u) := {u* G X* : u* is a subgradient of L at u} 

is called the subdifferential of L at w. Thus dL is a multivalued mapping defined on X and 
dL{u) = if L(u) = ±00. 
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Theorem 3.8. Let L : X \- > (—00, +00] be a proper convex and lower semi- continuous 
functional on the real Banach space X , then the subdifferential dL : X 1— > 2 X * is maximal 
monotone. 

Proof. See the fundamental paper by R. T. Rockafellar ( |14] ) . □ 

Theorem 3.9. Let X be a Hilbert space and A : D(A) C X n- 2 X . Then the following two 
statements are equivalent: 

1. A is a monotone operator and R(I + A) = X 

2. A is a maximal monotone operator. 

Proof. 

1. => 2. 

We only need to show that: if for any v G A (it), we have (v— v , u— Uq) > 0, then v G A{uq). 
Since R(I + A) = X, we can find U\ G X such that u + v = U\ + V\ for some V\ G A(ui). 
So (vi — Vq,ui — Mo) = —\\ui — u \\x > 0. Consequently, we have uo = u±, so vq G A{uq). 

2. =>- 1. See the fundamental paper by G. Minty([12]). □ 

Theorem 3.10. The operator A defined above is a maximal monotone operator on L 2 (fi). 
To prove this theorem, we need to introduce a p-Laplace type operator A p defined on L 2 (0): 

(u,v) G A p if and only if u G W 1 (L P (Q)) H L 2 {Q),v G L 2 (0) and 
J Q v(p dx = f n \Vu\ p - 2 Vu ■ V0 dx for any G H^(Z, p (n)) n L 2 (fi). 

From the definition of A p , we can see that v = A p {u) = — div(|Vit| p_2 Vtt) in T>'(Q) for 
it G 

Lemma 3.11. A p is a monotone operator and R(I + A p ) = L 2 (f2) ; i.e. A p is maximal 
monotone on L 2 (fl). 

Proof. 

(i) A p is monotone: let (ui,Vi), (ii 2 , v 2 ) G v4 p , then 

(«! - f 2, Mi - 1i 2 ) = / + t> 2 1i 2 - ViU 2 - f 2 «1 dx 

Jn 

= / |Vmi| p ~ 2 Vmi • Vui + |Vw 2 | p_2 Vw 2 • Vw 2 
Jn 

-|Vmi| p_2 Vmi ■ Vu 2 - |Vu 2 | p_2 Vu 2 ■ V«i dx 

> [ |Vwi| p + |Vw 2 | p - |Vmi| p_1 |Vm 2 | - |Vu 2 r _1 |VMi| dx 
Jn 

= [ (iV^r 1 - iVuar^dVuil - \Vu 2 \) dx 
Jn 

> 0. 
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(ii) R(I+A P ) = L 2 (fl): for any v G L 2 (fl), we need to find u G D(A P ) such that u+A p {u) 3 v. 
That is to say: given v G L 2 (fl), we need to find u G W 1 (L p (Q)) D L 2 (0) such that 

f(v-u)<f>dx= / |Vw| p - 2 Vw- V0 dx (3.6) 

for any <fi G W /1 (L p (fi))nL 2 (fi). It is actually a weak solution of the Euler-Lagrange equation 
for the following minimization problem: 

min T(u) (3.7) 

uew 1 (L p (n))nL 2 (fi) 

where T(u) := \\\v — u||i„ + -||Vm||^ . We will first prove there exists a minimizer for 
(13. 7p . Then we will show that such a minimizer is a solution for ( 13. 6p . Select a minimizing 
sequence {«fc}^ =1 for (13. 7p . Without loss of generality we can assume J n (v — Uk) dx = 
for every u k , otherwise if p- J Q (f — u fc ) dx = c 7^ 0, then \\v — + c)||| 2 < \\v — Wfe||| 2 - 

Let m = inf ueW i( ip (n))ni 2 (n)^(«) < 00. Since T(w fc ) ->■ m and || Vw fc |U P < pT{u k ), we have 
sup fc ||Vwfc||i < 00. Since f n Uk — u\ dx — 0, applying Poincare inequality, we have: 



||«/c||l p < \\Uk - Ux\\ Lp + < C\\Vu k - Vui\\l p + \\ui\\l p 

< C||VM fc || Lp +C||VMi|| Lp + HwiH^. 

So {uk}^i is bounded in W 1 (L P ('D,)). In addition, since ||w/e||z 2 — \f2T(uk) + |M|l 2 > we also 
have {uk}f? =1 bounded in L 2 (Q). Consequently there exists a subsequence {ix*-}^ C {uk}^ =1 
and a function u G W 1 (L p (f2)) such that u\ tj — ^ w weakly in W /1 (L p (f2)) and in L 2 (Q). So 
T(u) < liminfj_ i . 00 T(Mfc j .). By the fact that {uk^JLi is a minimizing sequence, we have 
T(u) < m. But from the definition of m, m < T(u). Consequently u is indeed a minimizer. 
Now we will show that u is a solution for (I3.6p . 

5T(u;<p) = Um-{T(u + e4>) -T(u)} 

e-S-0 e 

= Jim - (u + efiWl -\\v- u\\l 2 ) + i(||V(« + e0)f Lp - ||Vu||y } 

= lim - / (u - v)e(j) + \Vu\ p ~ 2 Vu ■ (eV0) + 0(e 2 ) 
(u - w)0 + I Vm| p_2 Vm • V0 dx. 

The necessary condition for u to be a minimizer of T(w) is ST(u; 0) = 0, which is (13.61) . □ 

Roughly speaking, we want to see u A(u) = lim^i A p {u) n . It is easy to show that A is 
monotone, to prove it is a maximal monotone operator on L 2 (f2), we also need the range 
condition: R(I + A) = L 2 (fi). We need the following two lemmas: 

Lemma 3.12. A is a monotone operator on L 2 (f2). 
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Proof. Let (ui,Vi), (u 2 ,v 2 ) G A. We have: 

(V 1 -V 2 ,U 1 -U 2 ) = + (f 2 ,W2) - <Vl,«2) - (V2jUl) 

f !«! + t> 2 M2 — ^1^2 — f 2^1 diC 



/ (zi, Dui) + (z 2 , -DM2) - (zi, -DM2) - (z 2 , Dux) 
Jn 

\ui\bv + Mbv - / (zi,Dtt 2 )- / (z 2l Du 1 ) 



Since -DM2) + (^2, Lh/i) < |«2|bv + I m i|bv ; we get ( v i ~ v 2, u i — u 2 ) > 0, which means 

A is monotone. □ 

Lemma 3.13. R(I + A) = L 2 (ft) 

Proof. We only need to show that: for any v G L 2 (Q), there exists u G BV(fi) fl L 2 (il) such 
that (u, v — u) G A. 

By lemma [3T2l we know that: given v G L 2 (Q), for any p > 1, there is w p G D 
L 2 (f2) such that (w p ,t> — u p ) G A p . Hence we have 



(v — u p )4> dx — |Vu p | p 2 Vu p ■ V(j) dx 
Jn ' Jn 



(3. 



for every G W^L^fi)) fl La (ft). 

Since A p is maximal monotone, by lemma ESI and noticing G A p (0), we have 

IKI| La = + A p )- l v\\ L2 < \\v\\ L2 . (3.9) 

Take <fi = u p and combine (I3.9p . we get the estimate: 

/ \Vu p \ p dx= / (v — Up)u p dx < \\v — u p \\l 2 \\u p \\l 2 < 2||u||| 2 — M, (3.10) 
Jn Jn 

for any p > 1 . 

By using Holder inequality we also have: 

/ |Vitp| dx<( ldx)«( \Vu p \ p dx)p = |fi|«Afp < max{|ft|M, |0|,M, 1} = M . 
in Jn 

Hence {w p } p >i is bounded in H /1 (L 1 (ft)) and we may extract a subsequence such that u p 
converges in Li(fi) and almost everywhere to some u G Li(Q) as p — > 1+. 
By Fatou's lemma, we have 

u 2 dx < liminf / u 2 dx < IM' 2 
n p-^+ Jn 



By lemma 13.11 we get 



Du\ < liminf / \Vu p \ dx < M , 
n ' p^ 1+ Jn 
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so we have that u G BV(f2) fl L 2 (fl). Since {u p } p> i is bounded in L 2 (fl), without loss of 
generality, we can assume u p — ^ u weakly in L 2 (£l) as p — > 1+. 

Now we want to show that {| Vu p \ p ~ 2 Vu p } p> i is weakly relatively compact in L2(f2;IR d ). 
First we need to obtain the following two estimates: 

( IV^- 1 dx<{[ \Vu p \ p dx)^^ < M^lttl? < M 
Jn Jn 

and for any measurable subset E C Q, 

| / \Vu p \ p - 2 Vu p dx\ < I \Vu p \ p ~ l dx < M^\E\v < max{M, 1}\E\% , for \E\ < 1 and 1 < p < 2. 

J E J E 

Hence {| Vu p \ p ~ 2 Vu p } p> i is bounded and equiintegrable in and consequently 

weakly relatively compact in L^fi;!^). Thus without loss of generality, we can assume: 
there exists z G Li(Q; R d ) such that, 

\Wu p \ p ~ 2 Wu p -± z as p ->■ 1+, weakly in Li(£l;R d ). 

Take <p G C^{Vt) in (EHD and let p -»■ 1+, we obtain 

(v — u)4> = z ■ V0, 
in in 

which means v — u = — div(^) in P'(fi). 

Now we need to prove ||z|| ioo < 1. For any k > 0, let B pk := {x G Q : |Vm p (x)| > A;}. By 
(13.101) . we have: 

< f \Vu p \ p dx< [ \Vu p \ p dx<M 



p,k 

Hence \B Pjk \ < § < max{^, §} for every 1 < p < 2 and k > 0. 

As above, there is some G such that | Vu p \ p ~ 2 Vu p Xb k ^ 9k weakly in 

as p -> 1+. 

Now for any G £oo(^; with H^Hl^ < 1, we can prove that 



/ \Vu p \ p - 2 Vu p X Bpk -4>dx\< [ iVMpI^" 1 dx< \ I \Vu p 



k 



Since 



/ 9k-<P dx\<\ / \Vu p \ p 2 Vu p Xb P:1! ■ <t> dx\ + | / (|Vm p | p 2 V« p XB p . fc - gk) ' <f> dx\. 
Jn Jn Jn 



let p — >■ 1+, we get: | f n g k • <f> dx\ < ¥■. So f n \g k \ dx < ¥■ for every k > 0. Hence g k — > i 
and a.e.. 
Since we have: 

| |Vu p | p - 2 VupXn\B p ,J < F" 1 < max{fc, fc 2 } for 1< p < 2, 



in 
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there exists fk G L tX) (Q;K. d ) and we can assume |Vw p | p 2 Vu p Xn\B k — ^ fk weakly in 
L 00 (f2;K c! ) when p -> 1+. Thus \\f k \\ Loo < liminf p ^ 1+ || |V« p | I '- 2 V« I ,Xn\B J ,.JUoc < L 

Since for any k > 0, we can write z = fk + gk, then we have z — fk = gk a.e.. By the 
fact H/fclUoc < 1, we get ||*|| Loo < 1. 
To prove (it, v — u) G A, we also need to show 

/ (0 - u)(v -u)dx= [ (z, D<f>) - \u\bv for any G W^L^tt)) n L 2 (fi). 

For any G W^L^fi)) nL 2 (fi), let n G C°°(fi) be such that <j) n in W^l^fi)) nLa(fi) 
as n — 7- oo. Using; </>„ Up clS cl test function in (13. 8p . we get 

(v- u p )((j) n - u p ) dx = |Vw p | p_2 Vw p ■ V(0„ - Up) dx. 
Jn Jn 

Hence, 

/ (v — Up) (0„ — Up) dx + / |Vw p | p dx = / |Vm p | p_2 Vm p • V0„ c?x. (3.11) 
Jn Jn Jn 

Since w p — u weakly in L 2 (fl) as p — > 1+, we have \\u\\l 2 < liminf p _s.i + ||w p ||l 2 , f n (v — 
u p )(pn dx — >■ Jq{v — u)(j) n dx and J n t>u p da; — >■ J n to da;. 

And since we have ||Vw p || Ll < || Vii p ||L p |f2| 1_ p and |m|bv < liminf p _ ill+ ||Vm p || Li , then 

|m|bv < liminf || Vup\\ L \VL\ l ~p < (liminf ||VwJ| l )( lim = liminf ||Vu p ||r ■ 

p— ►!+ p— ►!+ p— >l+ p— >1+ 



So 



|w|bv < (liminf || Vu p \\ L ) (lim sup || Vu p || L 

< (liminf ||VuJ|? )( lim M 1 ?) = liminf IIVuJI? . 

Combining (13. lip , we get: 

I (v — u)((p n — u) dx + \u\bv < z • V0„ dx. 
Jn Jn 

Letting n — > oo, we get: 

/ (v — u)(0 — u) dx + |u|bv < / z ■ (ix. 
Jn Jn 

By Theorem 13.41 for any <fi G BV(f2) fl L 2 (fl), we can find a sequence {0 n } C W /1 (L 1 (f2)) D 
L 2 (fi) such that: 

</>„->•</> in L 2 (Q) and / (z,D(f> n )^ / (z, £></>). 

in in 

So we can claim that 

(v — m)(0 — u) dx + \u\bv < / i z iD(f)) 
n in 
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holds for any G BV(f2) fl L 2 (fi). Consequently, we have: 

(v — u){4> — u) dx + |m|bv — / { z i D(p). 
n Jn 

So the theorem is proved. □ 
Recall' 

j( \ - ! Mbv u G BV(fi) n L 2 (tt) 

W : ~ { +00 ue L 2 (fi)\(BV(fi) n L 2 (ft)) 

J is a proper convex functional defined on L 2 (fi). Since (L 2 (fi))* = L 2 (fi), we have d J C 
L 2 (f2) x L 2 (fi). Now we can present the fundamental theorem in this subsection: 

Theorem 3.14. d.J = A. 

Proof. Since the functional J{u) is proper convex and lower semi-continuous, by Theorem 
13. 8[ we know dJ is maximal monotone. It is easy to see A C dJ . As we proved above, A is 
maximal monotone. So dJ = A. 

□ 

Remark. The results in this section also holds for Q = M. d . 

3.3 Characterization of the Pair (L 2 (Q), BV(Q)) 

From the definition of sub differential, we get the following Minimum Principle: 

Theorem 3.15. Let L : X \-t (—00, +00] be a proper functional on the real Banach space 
X . Then u* is a minimizer for L(u) if and only if G dL(u*). 

Since dT(u) = tdJ(u) — (f — u), the necessary and sufficient condition for u t to be a 
minimizer of problem (13. 2 p is: 



tdJ{u t ) - (f - u t ) 3 0. (3.12) 
Since we have shown in Theorem 13.141 that A = dJ, (13. 121) can be rewritten as: 

tA{u t ) - (/ - ut) 3 0. 
Theorem 3.16. The following assertions are equivalent: 

1. u t is a minimizer for problem A3.2\) . 

2. Ut G BV(fi) 2 and there exists z G X(Q) 2 with WzWl^ < 1 , f — u t = —tdiv(z) such that 



div(z)((p — Ut) dx — (z, D(j)) — \u t \ 
a Jn 



BV 



for any (f> G BV(fi) 2 . 
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3. Ut G BV(fi) 2 and there exists z G X(0) 2 with H^Hl^ < 1, f — Ut = —tdiv(z) such that 

f n (z,Du t ) = \ut\sv 
[z,u}=0 

Proof. 

2. =>■ 1. Since H^Hl^ < 1, we have f n (z,D<fi) < |0|bv- Then: 

-(/ - u t )((j) - u t ) dx < |0|bv - Nbv- 

Since 

^((/-^) 2 -(/-0) 2 ) = \U-^){<P-Ut) 

= -(/ - «t)(0 - Ut) - 7^(0 - U t ) 2 

< -rif ~ u t){(f> ~ u t ), 

we have ^||/ - mJ| 2 - ^||/ - 0||| 2 < |0| B v - KIbv for any <p G BV(1T) 2 . This tells us that 
Ut is a minimizer for problem (13. 2p . 

3. =>- 2. By Green's formula and the boundary condition [z, v\ = 0, we have 



div (z) (4> — u t ) dx = / (z,D(f>)— / (z,Du t ). 

Jn Jn 

Plug in J n (z, Du t ) = \ih\bv- We get - f Q div(z)(</> - «t) dx = j Q (z, D<j>) - |w 4 | B v- 

1. =>- 3. Since u t is a minimizer for problem (I3.2p . we have tA(u t ) 3 (/ — Mt). So there exists 

z G X(fi) 2 with ||-2||loc < 1, = f (/ — ^t) = — div(z) such that: 

J Q t>t0 = -D0) for any G BV(fi) 2 
J n (z,Du t ) = \ih\bv 

So (z, = Jo 0div(z) da; + J n (z, D<fi) = for any G BV($1) 2 . Hence [z, v\ — 0. □ 

Remark. From Theorem 13. 161 we can see that the Neumann boundary condition is a natural 
assumption for problem (13. 2p and f Q f — u t dx = J n — tdiv(z) dx = — J du t[z, v\ dH d ~ x = 
will automatically hold. 

In order to go further, we denote the homogeneous part of L P (Q) by L^(Q) := {v G 
L p (Q) : J Q v dx = 0} and define X(Q) p := {z G X(Q) p : [z,u] = 0}. As a rather deep 
result of [1], Bourgain and Brezis prove that for every v G L^(Q) with p > d, there exists 
2 G X(0) p such that v = div(z). This result is obviously not true for p < d and so we 
introduce the following norm for v G L^(Q) with 1 < p < oo: 

\\v\\ Y := inf {liminf \\z k \\ Loa : lim ||div(z fe ) - v\\ L =0, z fc G X(fi) p } 

and the corresponding normed vector space: 

y(fi) p := {v G LP(fi) : ||v|| yp < oo}. 
Then we have the following characterization for Y(Q) P : 
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Theorem 3.17. For every v G Y(Q) P , there exists z G X(Q) P such that 

v = div{z), \\z\\ Loo = \\v\\ Yp . 

In addition, the unit ball U p := {v G L^Q) : \\v\\ Yp < 1} is closed in L p norm topology. 

Proof. From the definition of the || • \\ Y norm, we can find a sequence {z k } C X(f2) p such 
that: 

lim \\div(z k ) - v\\ L = 0, lim \\z k \\ Loo = \\v\\ Yp . 

k— >oo fc— >oo 

Hence {H^Hl^} is bounded, so up to an extraction, we can find z G £<x>(^; such that z k 
converges to z weak-* in Loo(f2; R d ). Then for every G C°°(f2), 

v(j) dx = lim / div(zfc)0 dx = lim — / z k ■ V0 dx 

Jn k ^°°Jn k ^°° Jn 

= - I z • V0 dx = div(z)0 dx- [z, v]<i> dH' 
Jn Jn Jan 



d-l 



Choose G Cq°(Q), we get div(z) = t> G L p (f2) in the sense of distribution. So for any 
G C°°(n), we have: 

/ [z, u](j> dH^ 1 = 0. 
Jan 

Hence [z, u] — id d-1 a.e. on cftl 

By weak-* lower semi-continuity of || • || ioo , we get: 

\\z\\ t < lim \\zu\\ j = \\v\\y ■ 

II ll-^oo — , II K-ll-^oo II ll-tp 

By the definition of the || • \\ Yp norm, we have H^Hl^ > IMIv So H^Hl^ = ||f \\y p - 

Now let {v n } be a sequence in U p such that v n — > v for some v G Lq(J1), we want to show 
v G U p . Since w n = div(z n ) with z n G X(fi) p and \\v n \\ Yp = H^nlU^, we have H^Hl^ < 1- 
Thus we can find z G Loo(^; such that, up to an extraction, z n ^ z weak-* in Loo(^; 
So H^IIl^ < liminf^oo H^Hl^ < 1- For any G C°°(Q), we have: 

/ v(f) dx = lim / v n (f) dx = lim — z n ■ V0 dx = — z ■ V0 dx. 
in n ^°°Jn n ^°° Jn Jn 

Since — f Q z-V(j) dx = f n div(z)(j) dx — J dn [z,u](f) dH d ~ l , pick G C£°(J1), we get t> = div(z). 
Consequently, [z, z/] = 0. So f G £/ p . □ 

Lemma 3.18. For {v n }™ =l C and v G itfiJ/i: 

sup ||wn||L p < oo and lim \\v n — v\\ Yp — 0, 

we have v n — 1 t> weakly in L p (Q). 
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Proof. Since {v n }^ =1 is bounded in L P (Q), we can extract a subsequence {f refc }fc^i such that: 

v nk — ^ v weakly in L p (fl) (3.13) 

for some v G Lg(f2). Since lim^oo \\v n — v\\y = 0, without loss of generality, we can assume 
sup neN \\v n — v\\y p < oo. By Theorem 13 .17[ we can find {z n } C X(Q) p such that 

v n - v = div(z n ) and \\z n \\ Loo = \\v n - v\\ Yp . 

Since [z n , i>\ = 0, combining Theorem 13. 5[ we have: 

/ dw(z nk )(f) dx = - z nk • V0 dx + [z nk ,u]tr((j)) dH 
Jn Jn Jan 



d-l 



= - / Z nk -V(j) dx 

Jn 

for all (f) G fl L q (Q). Since lim^oo H-SnJUoo = 0, we have: 



lim / (f nfc — v)(j) dx = — lim / z rik ■ V0 <ix = 0. 

k^oo J n k^oo J Q 

This together with (13.131) shows v = v. Consequently, every subsequence of {v n }^ =1 has in 
turn a weakly convergent subsequence with limit v. This implies that v n — ^ v weakly in 
Lp(ft). □ 

With these preliminaries in hand, I am able to give two characterizations of the minimiz- 
ing pair (u t ,Vt). The first is the following. 

Theorem 3.19. Let (u t ,v t ) be the minimizing pair of problem 113. 1\) and f := ^ f n f dx. 

We have: 

(i) \\f-f\\y 2 <t^u t = f. 

(ii) ||/ - f\\y 2 > t O \\vt\\ Ya = t and J Q v t u t dx = t\u t \ BY . 

Proof. u t is a minimizer for problem (13.21) is equivalent to say (/ — u t ) G tA{u t ). So there 
exists z G X{Q) 2 with H^Hl^ < 1, f — u t — —tdW(z) in T>'{0) such that f Q (z,Du t ) = Im^bv- 
Since f n (z,Du t ) < II^IIl^ |^|bv < I m <|bv, the equality holds when (^Ibv — or H-zHl^ = 1- 

1. When I^Ibv = 0> we have u = constant. To minimize ||/ — Wt||i, 2 , we get u t = 

|^| J n f dx. SO ||/ - f\\y 2 = || - tdiv(s)||y a < t\\z\\ Loo < t. 

2. When |m(|bv > 0, we have H-zHl^ = 1- So we have 

/ v t ut dx = —t / div(z)ut dx = t (z,Du t ) = t\u t \BV- 
Jn Jn Jn 

And we claim ||t>t||y 2 = ^H-^Hioc = t. Otherwise, there exists z G X(Q) 2 with v t = 
— tdiv(z) and PHl,^ < 1, then J n v t u t dx = t f n (z,Du t ) < t\u t \Bv, which contradicts 
our previous statement. 
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□ 

Remark. The above theorem is a generalization of Yves Meyer's result in [TTJ which was 
proved for the special case Q = M. 2 by using techniques from harmonic analysis. 

Antonin Chambolle [6] has introduced another type of characterization for a finite di- 
mensional minimization problem related to (13 . 1 j) . I have generalized his result to the case 
(L 2 ,BV) as a consequence of theorem 13.191 

Theorem 3.20. Given f G L 2 (Q), the minimizing pair for problem ^3. 1\) is (ut,Vt) = 
(/ ~~ TTtt/jjC/Oj 7r %(/))j where ittu 2 (f) the L 2 projection of f onto the set tU 2 . 

Proof. 

1. When ||/ - /||y 2 < t, we have u t = n tU2 (f) — f — f E tU 2 . 

2. When ||/ — f\\y 2 > t, by Theorem 13.191 we have ||ft||y 2 = t and J n v t u t dx = t\u t \BV- For 
any w G tU 2 , there exists z G X(0) 2 such that w = —div(z) and U-zHl^, = HHta — Thus 

10 (/ ~ u t) dx — wilt dx = l — div(^)u i dx = (z,Du t ). 
n Jn Jn Jn 

" " ^ " |BV- 



Thus we have: 



/ w(f - v t ) dx < II^IUooI^Ibv < t\u t \ 
Jn 

(w — v t )(f — v t ) dx < t\u t \BV — t\u t \BV = for any w G tU 2 . 

Jn 

So we have v t = rr tU2 (f) and u t = f - nu 3 {f)- □ 

Since we have characterized the minimizing pair (ut,Vt), we can now give an alternative 
expression for the K-functional: 

K(f,t) := inf {i|H|| 2+f |u| BV }. 

f=u+v I 



Theorem 3.21. 



K(f,t)= j W/)/^-^IIW/)!IL- 



Proof. From Theorem 13.191 we know that (u t ,v t ) — (f — 7Ttu 2 (f),^tu 2 (f)) an d fn v t u t dx = 

t\u t \BV- So 

K(f,t) = \htu 2 (f)\\l 2 + ju - 7rtu 2 (f))nu 2 (f) dx = J ntuMf dx - ^||W/)llL- 

□ 

By using Theorem 13.161 we can calculate minimizers explicitly for some simple cases. 
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Example 3.22. Let Q = B(0,R) = {x e R d : \x\ < R} and / = Xb(o,v) with < r < R. 
When < t(\dB(0,r)\/\B(0,r)\ + \dB(0,r)\/{\B(0, R)\ - |£(0,r)|)) < 1, we have: 

d d ■ r d_1 

U t = (1 ~ t-)X B {0,r) + t Rd _ rd X B (0,R)\B{0,r) 

\dB(0,r)\ \dB{0,r)\ 
~ (1 " 1 \B(0,r)\ )Xmr) + V(0,i?)|-|fi(0,r)| XB(0 ^ (0 ' r) 

and 

d d ■ r d_1 

Vt = t-X B (0,r) ~ t Rd _ rd XB(0,R)\B(0,r) 

. \dB(0,r)\ \dB(0,r)\ 
~ 1 \B(0,r)\ XB(0 ' r)_ V(0,i?)|-|5(0,r)| XB(0 ' RAB{0 ' r) - 

When t(\dB(0,r)\/\B(0,r)\ + \dB(0,r)\/(\B(0, R)\ - \B(0,r)\)) > 1, we have: 

Ut = Rd XB(0 ' H) 

W,r)\ 
\B(0,R)\ 

and 

Vt = X B (o, r ) — -j~[Xb(o,r) 

~ Xmr) ~ \B(0,R)\ XmR) - 

Proof. We look for the minimizer u t with the form u t = aX B (o, r ) + bXB(o,R)\B(o,r)- Then 
tdiv(z) = -(/ - u t ) = (a - l)X B (o, r ) + &Xs(o,R)\B(o,r) ■ We take z = ^jp-x for x E -6(0, r), 
then Ut = f + tdiv(z) = a for x G -8(0, r). To construct z in -5(0, R)\B(0, r), we will look for 
with the form z = p(\x\) o- Since H^Hl^ < 1, we need p(r) = — 1, this tells us 



Xb(o,r) 



z 



_ _i 



So a = 1 — — . To make [z, v\ = 0, we require p(R) = 0. Since 

div(z) = Vp(|x|) • ^ + p(|x|)div(|| ) = p'(|x|) + A>(l*l)^p 

we must have: 

t(p'(s) + p(s)*=±) = b for r < s < R 
p(r) = -1 
p(R) = 

Solve this ODE, we get: p(s) = -^fz^s 1 ^ + -j^s and b = ffarz^- In addition, p'(s) = 
((d - l)(f ) d + > 0, so -1 < p(s) < for r < s < R. 

Thus, U t = (1 -i;)Xfi(0,r) + t ^_ r d X B (0,R)\B(0,r) and V t = tfXB(0,r) - t Rd_ r d X B (0,R)\B(0,r) ■ 
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To show u t is a minimizer for (13. 2p . we only need to check whether f n (z, Du t ) = \ut\sv- 
By Green's formula, we have: 



/ (z, Dut) = — I div(z)ut dx — — I v t u t dx 
Jn Jn t Jq 



1 r f ,-> d^ d 1 f d-r^ 1 . d-r d - 1 . , , 

= ?{ (! - dx + / fc(-fe) 

1 JB(0,r) r r J B(0,R)\B(0,r) n 7 n r 

d d d- r d ~ x d ■ r^ 1 

= ^-t-)-\mr)\-t Rd _ rdRd _ rd (\B(0,R)\ - |S(0,r)|) 

= (l-t±-t±^)H<-\dB(0,r)) 

= \ut\BV- 

The above equality makes sense only when 1 — 1(- + j^r—fx) > 0, which means 
< *(|9S(0,r)|/|S(0,r)| + \dB(0,r)\/(\B(0, R)\ - \B(0,r)\)) < 1. 



□ 



Example 3.23. Let f(x) = x for x G [0, 1]. 

For < t < |, we have: 

( V2i < x < V2t 
u t = I x V2t <x<\- \plt 

{ l-y/2t I - V2t <x <l 

For t > |, we have: u t = \ for x G [0, 1]. 

Proof. Since te' = — v t and H^Hx^ — 1 3 we look for z with the following structure: 

-j- t (x-h) 2 + l 0<x<h 
1 h<x<l-h 
~Y t (x- (1 - /i)) 2 + 1 l-/i<x<l 

We need to choose h such that z(0) = z(l) = 0(Neumann Boundary Condition). So we have 
1 — ^ = 0. Hence h = y/2t. But we also require h < ~, so \/2i < |. 
Thus for t < i we have: 



-j- t (x - V2t) 2 + 1 < x < v 7 ^ 

1 v 7 ^ < x < 1 - V2i 

l-y/2t)) 2 + l 1- V2i<x <1 



2t ( X 



Then we have: 



v 7 ^ < x < v 7 ^ 
//,. = <; x v 7 ^ < x < 1 - V2t 

1 - V2l 1 - V2t < x < 1 
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For t > |, we have z = — 4(x — \) 2 + 1 for x G [0, 1]. Then we have: ut = \ for x G [0, 1]. To 
show u t is a minimizer, we only need to check f Q (z, Du t ) = \u t \BV- 

1 nl-h 

(z, Du t ) — I z dx = 1 — 2h — I^Ibv- 
Jh 

□ 



4 Multiscale Decompositions 

The solution of minimization problems like ( 13. ip leads to multiscale decompositions of a 
general function /. In the case we have been considering, each / G L 2 (fl) is decomposed as 
=0 Wk where each w k is viewed as providing the detail of / at some scale. Currently, 
there are several ways to achieve this goal. The most common of these is to use a standard 
telescoping decomposition where w k := u tk — "Ut i ._ 1 and t k = 2~ k . Other approaches to obtain 
multiscale decompositions were given by Eitan Tadmor et al.'s work (see [16]) and Stanely 
Osher et al.'s work (see [13J). 

Since ||ft||y 2 for problem (13.11) depends on the parameter t, this gives us a way of decom- 
posing a given function / G L 2 (Q) into different components based on the size of the || ■ \\y 2 
norm of each component. This approach falls into a category of methods (called Inverse 
Scale Space Methods) that were introduced by Groetsch and Scherzer in [9]. It centers on 
using the above || • ||y p norm to measure the oscillation of v t in a cetain sense. In our language, 
the choice of components takes the following form: 

u k+ i := argmin {- 1|/ - u\\l 2 + t k J{u, u k )}, (4.1) 

where |||/ — w||| 2 is the L2-norm fit-to-data term and J(u,u k ) is a regularization term. 
Typically we initialize uq = or uq = pr f n f dx and we require that {u k } satisfies the 
inverse fidelity property: 

lim ||/ - u k \\ L2 0. 

If, as a special case, we consider BV minimization and choose J(u, u k ) as the Bregman 
distance defined by 

D(u,u k ) := |u|bv - KIbv - / s(u — u k ) dx, 

Jn 

where s G d(\u k \-Bv), then ( 14. ip becomes the method introduced by Osher et.al. in [13] . This 
method has many promising properties for image denoising which were proved in [13"] . 

The interpretation of the Bregman distance for image processing given above is somewhat 
ambiguous. Another possibility is to simply take J(u,u k ) := \u — u k \-Q\- Roughly speaking, 
\u k+ i — u k \-Qv measures the similarity between two images u k and u k+ ±. For any choice t k > 0, 
u k+ i contains more detail than u k and is closer to /. We choose a sequence t > t± > t 2 > ... 
with lim^oc t n = 0. One sees that 

u k+1 := argmin {-||/ - u\\l 2 +t k \u - u k \ B y}, (4.2) 
ueBV(fi)nL 2 (tt) 1 
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with wo = 0. If we take w k+ i '■= u k+ i — u k , wo := «o and v k '■= f — Uk, then (14. 2 p can be 
viewed as: ^ 

(tUfc+ijUfc+i) = argmin {-||v||£ 2 +t k \w\ B v}, 

with Vq = f. Thus Uk+\ = Yl^=i w n will be a minimizer for (14. 2ft . This is the hierarchical 
(L2,BV) decomposition method introduced by Tadmor et.al. in [16]. 

Theorem 4.1. Le£ / G ^2(^) and {uk}^^ be defined as in Then we have: 

1. ||/-u*+i|Ua < \\f - u kh 2 - 

2- ||/-Mfc+i||y a <t h ->0. 

3. Let v n := f - u n , then Y2=o{ 2t k\ u k+\ ~ W/c|bv + IK+i - Wfc||i 2 } = ||/||| 2 - IK+i||| 2 - 

Proof. 

1. Since 

11 1 

2 11/ ~ Mfe HL ^ _ Mfc+1 '^2 + tk \ Uk + l ~ Mfc l BV - 2 HZ ~ Mfc+1 Hi2' 

we have 

HZ-Mfc+ilk < ||/-«fe|U 2 - 

2. By Theorem 13. 161 

3. If \\v k - c k \\ Y2 < t k , where c k := p f n v k dx, then - u k ,v k+1 ) = (c k ,v k - c k ). 
Otherwise, we have: 

|K+i||y 2 = t k , (v jfc+i, u k+x - u k ) := / w fc+ i(M fc+1 - t^) = tfcl^fc+i - w&|bv- (4-3) 
Since - u k + Wfc+i = we get: 

IK Hz* = (( u fc+i-^fc)+^fc+i ; ( M fc+i- M fc)+^fe+i) = \\v k+ i\\l 2 + \\u k+1 -u k \\l 2 +2(v k+1 ,u k+1 -u k ). 

(4.4) 

Combining (14.31) and (I4.4p . we get: 

IKIlL ~ IK+iIlL = IK+i ~ u k\\\ 2 + 2t k \u k+1 - u k \ BY . (4.5) 
Sum (14. 5 P from k = to k = n, we get: 

n 

^{2t fc |M fc+ i - U k \ B V + \\u h +l - U k \\l 2 } = \\f\\l 2 - ||Vn+l|li a - ( 4 - 6 ) 
fc=0 

□ 

In addition, we have the following L 2 convergence result: 

Theorem 4.2. Let f G L 2 (Q) with Q bounded Lipschitz domain in M. d and t k = to ■ r k with 
< r < 1, then we have: 
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1. lim fc _ > . 00 ||/ - u k \\ L2 = 0. 

2 - EfcL { 2t fcK+l ^ U k\BV + \\u k+ l ~ U k \\l 2 } = \\f\\l 2 . 

Proof. 

1. By Theorem 14.11 we know that {||vn||L 2 } is a decreasing sequence. Hence, to prove 

IknlUa °> We OIll y need to SnOW ll^n+llUa -> 0. Note V 2n +1 = V n - YJkLn^k+l ~ u k) ■ 

Multiply V2n+i with itself, we get: 

2n 

\\v2n+l\\ 2 L 2 = -(Uan+i, J^(Ufc+l - Ufc)) + («2n+l,Un) == A + B. 

k=n 

By Theorem 14.11 we know ||f2n+i||y 2 < hn- So 

2n 2n 2n 

|^4| < hn\ ^(Uk+l - M/c)|BV < hn |«fe+l - «fc|BV < ^Jtfc|«fc+1 - Ufc| B V- 
fc=n fc=n fe=n 

From Theorem O, we know that J2k=o t k\ u k+i ~ w/cIbv < §||/||i 2 - Hence ELo^K+i ~ 
Uk\~Bv} is a Cauchy sequence. So |A| — > when n — > oo. 
Since v n = f - u n = f - Y2=ai u k+i ~ we have: 

n-l 



1 5 1 = | (f 2n +l, /) - 2n+l, ^ fc+i - Uk) \ 

k=0 

n— 1 

< |(^2n+l) /)! + *2n ~ M fc | B y 



k=0 
n-l 



t 

- fc=U 
^2n II n || 2 



< K^n+lj /)| + ^T" |M \\L 2 

Since lim n ^ ^ = 0, we have \B\ — ^ iff \(v 2n , f) \ ->• 0. 

Since m^+i is a minimizer for (I4.2p . there exists Zk+i G X(fi) 2 with || z k +i Wl^ < 1, ffc+i 
/ - u k +i = -t k div(z k+ i) = div(-t k z k+1 ) e L^(Q) such that 



f Q (z k+ i, D(uk+i - u k )) = \uk+i - u k 
[z k +x, v] = 



BV 



Since we have sup fcgN ||ffc||i 2 < \\f\\h 2 an d H^fc+ilta — ^k — >• 0, by Theorem 13.181 we know 
v n — ^ weakly in L 2 (tt). So we have \(v 2 n+i, f)\ — > 0. Consequently H^n+ilUa ~* 0- So 
lim^oo ||/ - Uk\\ L2 = 0. 

2. Recall from Theorem 14.21 we have: 



^2{2t k \u k+ i - U k \ BV + \\u k+1 - U k \\ 2 L2 } = \\f\\l 2 - HVn+lllia- 



fc=0 
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Let n — > oo and notice ||fn||z, 2 — 0> we have: 

oo 

^{2t fc |w fc+1 - M fc | BV + IK+i - u k \\l 2 } = \\f\\l 2 

k=0 

□ 

Remark. In [16] (Tadmor et al.), the same result was proved under the assumption / G 
BV(fi)(fi C R 2 ) or / G (L 2 ,BV)e with < 6 < 1. Here we have removed the smoothness 
assumption. 



JJu) :-- 



5 Decomposition for (X,Y) = (L pi W l (L T )) 

Now let's consider the following (L p (0), W /1 (L T (fi))) decomposition with 1/r := l/p + l/d: 

(u t , v t ) := argmin + *|| Vw|| Lt }, (5.1) 

u+v=f P 

where || Vu||x, T := (J n | Vu| r dx) 1 / 7 " and 1 < p < oo. In this section, Q := M. d or is a bounded 
Lipschitz domain in M. d . Thus by Sobolev embedding theorem, we have W 1 (L T (Q)) C L p (f2). 
Let 

\\Vu\\ Lt ueW\L r (ti)) 
+oo «GL p (0)\r(L r (0)) 

be a functional defined on L P (Q). Since (L p )* — L q for i + ^ = 1, by Riesz representation 
theorem, any functional s defined on L p (f2) can be represented as (s, v) := J Q st> dx for 
any t> G L P (Q). To study the pair (L p (Q), W 1 (L T (Q))), we first need to characterize the 
sub differential of the functional J T . For s G Lfj(Q) (i + - = 1), we can define the norm 

II • 1 1 Gt by 

H s IIg t := sup / sv dx 

Theorem 5.1. 

(i) For u 7^ constant, d.J T {u) := {s G Lq(Q) : f n su dx = ||Vm||l t and ||s||c? T = 1}- 

(ii) For u = constant, dJ T (u) := {s G L^(Q) : ||s||g t — !}■ 

Proof. Given a function u G H /1 (L r (fi)), we can define a functional s on M.u such that 
(s,cu) = c||Vm||l t for any c G M. By Hahn-Banach Theorem, we can extend the domain 
of the functional to W l (L T (Q)). Let's say functional s with s|r u = s and |(s, u)| < ||Vi>||l t 
for any v G W l (L T {VL)). Since s|r u = s, we have (s,u) = ||Vw||z T . Consequently, we have 
\\ s \\g t < 1 f° r s £ -^□(fi) (equality holds when Vu ^ 0). 
Then if s G Lg(fi), we have: 



(s,v — u)= / s(f —u) dx < \\s\\ Gt \\S7v\\ Lt - ||V«||i T < J T (v) - J T 



u 
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for any v G W x {L r {Q)). So s G d.J T {u). 

Conversely, if s G dJ T (u), then J T (v) — J T (u) > (s,v — u) for any v G W 1 (L T (Q)). In 
addition, we have s G L g ($7). By taking v = Xu, we get 

(l-A)((s,u)-||Vii|k) >0. 

By successively taking A > 1 and A < 1, we deduce that (s,u) = ||Vw||l t . Therefore, (s,v) < 
||Vv||l t for all v G W 1 (L T {VL)) and (s,u) = ||Vw||l t . This implies that ||s||g t < l(equality 
holds when Vu ^ 0). Since (s, c) < J r {u + c) — J T {u) = for any constant c, we have 
J n s dx = 0, which means s G L^(Q). □ 

Define the duality mapping $ p : L P (Q) \-t L q (fl) (- + - = 1) by: 

Z p (u) := \u\ p ~ 2 u. 

It is easy to check Z p (u) = d(-\\ ■ \\ P i p )(u). Then we have the following theorem for the 
minimizing pair (L P (Q), W 1 (L T (Q))). 

Theorem 5.2. Given f G L P (Q) and let (u t ,v t ) be the minimizing pair of problem ( 15. 1\) and 
c f := argmin cgR ||/ - c\\ Lp . We have: 

1- ||3 P (/-c/)||Gf T <t&u t = c f . 

2- \\3 P (f - Cf)\\ 0r > t <^> ||3 p (vt)||o T = t and J n Z p {v t )u t dx = t\\Vu t \\ LT . 

Proof. Since (3 p (f - cf, c) < - c f + c)\\ p Lp - \\f - c f \\ p Lp ) for any c G R and the right 
handside of the inequality is always nonnegative, we get (3 p (f — c/),c) = for any c G R, 
which means f n 3 p (f — cf dx = 0. So 3 P (f — cf G Lq(O). 

w t is a minimizer for problem (15.1 1) is equivalent to say Z p (f — u t ) G tdJ T {u t ). So Z P {vt) G 
L^(^). In addition, we have ||3pK)|| Gr < t and f Q Z p (v t )u t dx = t\\Vu t \\ LT . 

1. When ||3p(/ - c/)|| Gt < t, we have - cf G tdJ T {cf). So u t = c/. 

2. When \\3 p (f — Cf)\\a T > t, u t cannot be a constant. So ||3p(^t)llG T = t. 

□ 

6 Decomposition for (X,Y) = (tiAp) 

A special case that is important in analysis and in numerical methods is when X and Y are 
a pair of l p spaces. Such problems occur when we discretize the decomposition problems 
for Sobolev or Besov spaces and also when we develop numerical methods. In this chapter 
we shall study the minimizing pair for the case of X = £2 '■= ^(Z) and Y = £ p := f p (Z), 
1 < p < 00, i.e. the problem 

(x t , y t ) := argmin { ~\\y\\% + t\\x\\ ep }, (6.1) 

x+y=b ^ 

where b G £2- 
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Theorem 6.1. 

1. For x y^O, <9(||x||^) := {s e £ q : s ■ x = \\x\\t p and \\s\\ £q = 1}. 

2. Forx = 0,d(\\x\U p ):={se£ q : < 1}. 

Proof. Given a sequence x G £ p , we can define a functional s on Ma; such that (s, cx) = c\\x\\e 
for any cGl By Hahn-Banach Theorem, we can extend the domain of the functional to 
£ p . Let's say functional s with s\^ x = s and \{s,y)\ < \\y\\e p for any y G £ p . Since s\^ x = s, 
we have (s,x) = \\x\\e p . So ||s||^ < 1 and we have \\s\\e q = 1 for the case x ^ 0. Then we 
have 

s ■ (y - x) < \\s\\ eq \\y\\e p - \\x\\e p < \\y\\e p - \\x\\e p . 

Sosed(\\x\\ tp ). 

Conversely, if s G then \\y\\e p — \\x\\i p > s ■ (y — x). By taking y = Xx, we get 

(1- X)(s-x- \\x\\ £p ) > 0. 

By successively taking A > 1 and A < 1, we deduce that s-x — \\x\\i p . Therefore, s-y < \\y\\i p 
for all y G £ p and s-x = \\x\\e p . This implies that \\s\\e q < 1 (equality holds when x ^ 0). □ 

For the case p — 1, the minimizing pair (x t ,y t ) can be obtained by the "soft threshold- 
ing" procedure which is widely used for wavelet shrinkage in image processing (See [7]). 
That is to say: x\ = sign(6 J ) max{0, \b l \ — £}, where x t = {x\} and b = {&*}• It can be shown 
that the "soft thresholding" technique is a special case of the following characterization. 

Theorem 6.2. Given b G £2, let (xt,yt) be the minimizing pair of problem \6. We have: 

1. \\b\\ iq < t & x t = 0. 

2. \\b\\i q > t <^> \\y t \U q = t and y t • x t = t||^||^,. 

Proof. (x t , yt) is the minimizing pair of problem ( 16. lft is equivalent to say: b—Xt G £d(||a; t ||^J. 
Consequently, we have y t ■ x t = (b — x t ) ■ x t = t\\x t \\i p . 

1. When \\b\\ tq < t, we have b G td(\\ ■ \\e p ){0). So x t = 0. 

2. When \\b\\e > t, x t cannot be the zero element. So \\yt\\e q = \\b — x t \\e q = t. 

□ 

Let's define the convex set 

U q := {y G £ 2 : ||y||^ < 1} 

where q is the dual index to p (1/q + l/p — 1). Then the set U q is closed in £ 2 norm topology, 
so the £2 projection of a given sequence b G £2 onto the set U q is always well defined. 

Theorem 6.3. Given b G £2, the minimizing pair for problem \b\ 1\) is (x t ,yt) = (b — 
ir tUq (b) , ir tUq (b)) , where ^tu q ip) ^ s the £ 2 projection of b onto the settU q . 
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Proof. 

1. When \\b\\ iq < t, we have y t = Tc tUq (b) = b G tU q . 

2. When ||6||^ 8 > t, by Theorem 16.21 we have \\yt\\e q = t and y t ■ x t = t\\xt\\e p . For any 
z t G tU q , we have: 

(zt - y t ) • (b - y t ) = z t ■ x t - y t ■ x t < \\zt\\e 9 \\xt\\e p - t\\xt\U P < 0. 

So y t = iTtu q {b). Consequently, x t = b - n t u t (b). □ 

Remark. This result paved a road for studying (L 2 , Bp(L p )) type decompositions, where 
Bp(Lp) is a Besov space. Since the norms of L 2 and Besov spaces can be characterized 
by their wavelet coefficients, the (L 2 , Bp(L p )) type decompositions can be derived from the 
decompositions for sequence spaces. 

7 Numerical Implementation for (X, Y) = (L 2 (Q), BV(fi)) 

The problem of finding a minimizing pair (u t , v t ) almost always is solved numerically. Typi- 
cally, numerical methods are built through some discretization of the continuous problem. In 
this chapter, we will study the numerical implementation for (L 2 (fl), BV(fi)) decomposition. 
Without loss of generality, we can assume Q = [0, l] d . Let 

D n := {2- n ([k ll h + l]x [fc 2 , k 2 + 1] x • • • X [k d , k d +l}) : < k t < 2 n - 1, 1 < i < d} 

be dyadic cubes with length 2~ n . We replace / by a piecewise constant approximation of / 
which has the following form: 

/„ := ^2 a iXi- 

Thus, f n is a function in the linear space 

V n :— {/ G L 2 (fl) : / i s constant on /, I G -D n }. 

The spaces {V^}, n > 0, are nested, i.e., Vo C • • ■ C V n C C • • -. As our approxima- 
tion, we will take f n = P n (f), where P n (f) is the L 2 projection of / onto V n . The original 
problem (13. 2p is then approximated by the finite-dimensional problem 

u n := argmin {^||/ n - u\\ 2 L + *|u|bv}- (7-1) 

uev n * 

We can compute ||/ n — w n |U 2 and |m„|bv exactly by discrete norms. Problem ( 17. lft then 
becomes a ^x-minimization problem: 

x t := argmin {||6 - + t||Mx||^}, (7.2) 

a; 

where M is an m x n matrix with the property M(x + c) = Mx for any constant vector c. 
Define 

K n (f,t):= inf {l||/ n -„||2 2 +tH BV }. 

We have the following result about the convergence of the discrete solution to the con- 
tinuous solution: 
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Theorem 7.1. Let ut be the minimizer of problem 113.2]) . We have: 

1. u n — ^ ut weakly in L 2 . 

2. u n —> u t strongly in L p , where 1 < p < d/(d — 1). 

3. ]im^ 00 K n (f,t)=K(f,t). 

Proof. Since |||/ n - u n \\\ 2 < \\\f n - u n \\\ 2 + t\u n \ BV < |||/n||| 2 , we have \\u n \\ L2 < \\f n \\L 2 - 
Similarly, we can get \u n \ B y < ^||/n||| 2 - Since f n = P n (f) -»• / in L 2 , {||Mn||L 2 } and {|u n |Bv} 
are bounded sequences. So we can find a subsequence {u nk } such that: 

1. u nk — ^ u weakly in L 2 . 

2. u nk — > u strongly in L p for 1 < p < d/(d — 1). (by the compact embedding BV(O) e — >■ 

L p (n)) 

3- |w„Jbv — >" I^Ibv- 
for some u G BV(fi) fl L 2 (£l). Then we have: 

K(f,t) = ~\\f-u t f L2 +t\ut\ BV < ~\\f-uf L2 +t\u\ By < liminf K n „{f,t). (7.3) 

Given the minimizer w t G BV(f2) fl L 2 (fi), we can find a sequence {tiv} with w n G V„ such 
that: 

w n -)■ Mt in L 2 and |w„| B v ~> |«*|bv- 
Since K n (f,t) = |||/ n - w n ||| 2 + |u„| B v < |||/ n - w„||| 2 + Wbv, we have: 

limsvp K n (f,t) < limsup {-||/ n - Wn\\l 2 + Wbv} = K(f, t). (7.4) 

Combining (17.31) and ( 17.4ft . we know that u is a minimizer for problem ( 13.21) and 

lim sup = liminf K nh (f,t) = K(f,t). (7.5) 

n— >oo fe— >oo 

By the uniqueness of the minimizer, we have u = u t . Since we can extract a further subse- 
quence {u nk }, which satisfies (17. 5p and has the limit u t , from any subsequence of {u n }, we 
have: 

1. u n — ^ «t weakly in L 2 . 

2. w n — )■ -u t strongly in L p , where 1 < p < d/(d — 1). 

3. lim n ^ 00 K n (f,t) = K{f,t). 

□ 
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